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Abstract 

The second and third factorial moments of the multiplicity distributions of 
gluon and quark jets are calculated up to the next-to-next-to-next-to lead- 

> 

ing order in perturbative QCD, using the equations for generating functions. 

■ The results are confronted with experimental data. A general discussion on 

o ' 

■ high order corrections revealed by such an approach is given. Other possible 
r~| ■ corrections and approaches are discussed as well. 

a: 

^ . I. INTRODUCTION 

> : 

•rH , 

The detailed properties of multiplicity distributions of gluon and quark jets are stud- 
ied in experiment nowadays [IIH3]. Theoretical description of the data in the framework of 



perturbative QCD is rather successful qualitatively, often at the accuracy of 10-15% or bet- 
ter. Such an accuracy, obtained sometimes only in low-order approximations, is somewhat 
surprising because the QCD expansion parameter is rather large at present energies (about 
0.5), so higher order contributions should be estimated. It happens that different physical 
quantities are sensitive in a different way to higher order corrections. A thorough discus- 
sion of average multiplicities of gluon and quark jets is given in [4-8] with respect to their 
dependence on these corrections. Here we extend these studies to the widths and higher 
moments of multiplicity distributions of gluon and quark jets, calculating them in the next- 



1 



to-next-to-next-to leading order approximation. Such an accuracy becomes possible because 
the equations for the generating functions admit a perturbative expansion. Besides finding 
out the appropriate coefficients, we also discuss some general problems, revealed by such an 
analysis of average multiplicities, widths, and higher moments. 

II. WIDTHS OF MULTIPLICITY DISTRIBUTIONS IN PERTURBATIVE QCD 

Any moment of the parton multiplicity distribution of gluon and quark jets in QCD can 
be obtained from the equations for the generating functions 

G' G = f 1 dxKg(x)>y*[G G (y + \nx)G G (y + ln(l - x)) - G G (y)\ 
Jo 

+ n f f dxK%(x)>y*[G F (y + \nx)G F (y + ln(l - x)) - G G (y)}, (1) 

J 

G' F = f 1 dxK^(xh 2 [G G (y + \nx)G F (y + ln(l - x)) - G F (y)], (2) 

where Gi are the generating functions of the multiplicity distributions of gluon (i = G) 
and quark (i = F) jets, defined by 

oo oo q 

G t (y,z) = ^> + l) n ^ } = E^W^ W - (3) 

n=0 q=0 

In these expressions, (n*) = Y.^Lo n Pn^ is the average multiplicity, z is an auxiliary variable, 
y = ln(p©/(5o) is the evolution variable, p, 6 are the momentum and the opening angle of 
a jet, Qo=const , G'(y) = dG/dy, and rif is the number of active flavors. Moreover, 
2 2N c a s {y) 2tt / (3 1 \n2y \ 

117V. - 2n f l7N*-n f (5N c + 3C F ) 
Po = g , Pi g , (oj 

«5 is the running coupling strength, N c is the number of colours, and C F = (N 2 — 1)/2N C = 
4/3 in QCD. The argument of 7g in the integrals is chosen to be y+lnx(l — x), as determined 
by the transverse momentum of partons at the splitting vertex. The kernels of the equations 
are 



Kg(x) = - - (1 -x)[2-x{l -x)}, (6) 



K F G {x) = ^-y + (l-x)\ (7) 



K°(x) = ^ 



--l + X - 

x 2 



The normalized factorial moment of any rank q can be obtained by differentiation 



(8) 



F (i) = 1 d9Gi \ . (g) 

or, equivalently, by using the series and collecting the terms with equal powers of z on 
both sides of the equations (jl]), (@). 

The normalized second factorial moment F2 defines the width of the multiplicity distri- 
bution, and is related to its dispersion D 2 = (n 2 ) — (n) 2 by the formula 

D 2 = (F 2 - l)(n) 2 + (n) = K 2 (n) 2 + (n), (10) 

where K 2 is the second cumulant. 

The second factorial moments normalized to their own average multiplicities squared are 

jpG _ ( n G(n G - 1)) fF _ (n F (n F ~ 1)) nn 
2 — / \2 ' 2 — / \2 ' v / 

{n G y (n F y 

Let us write down their perturbative expansions up to 7g-terms as 

^ 2 G = ^(l-/i7o-/27o 2 -/37o 3 ), (12) 

Fi = (1 + ~ 0i7o - 027o - 037o 3 ), (13) 

where r = N C /C F determines the asymptotical value of the ratio of multiplicities in gluon 
and quark jets (see (|T8| ) below). It is equal to 9/4 in QCD and to 1 in SUSY QCD. Actually, 
the asymptotical values of F 2 and F 2 in front of the brackets in (|T2"D, fll3]) are found out 
from the equations below by equating the leading terms on both sides. However, we have 
inserted their explicit expressions directly here to simplify further notations. 



Using the Taylor series expansion of G's as proposed in 0, one can rewrite the Eqs. (|J), 
]) as 

\[lnG G ]" — Gq — 1 — 2G' gVi + G' G v 2 + 0.5Gv> 3 + 
7o 

(If ) Vl2 + 4jB7 o[ (Gg " 1)vi " G ^ 2 ] + 
/ , , / , ,, / , / 

lit rl Lt j? \ _ / upLrc \ / UfUp \ / Lr 



1 — V 4 + 2 — V 5 + — y 6 + — V13 



4N C \ Gq ) V Gg* / \ J V Cg , 



2B7„M(^-l k + 2 fgE) U|], (14) 



^[lnGV]" = G G - 1 - Gg v 7 - G G v 8 - 0.5G G v 9 - 
7o 

/ G G G F \ i G G G F \ ( G G G F \ 

2fi 7o 2 [(G G - l)v 7 + 2G G v 8 + (G' G + %^)v u ]. (15) 

The terms up to the third derivative of G are kept everywhere because each derivative gives 
rise to the factors 7 or 7^ as seen from their definition below in Eq. fll6|) , and we are 
interested in corrections up to 7^. We use B = j3 /8N c , B\ = /3i/4/3 N c . The corresponding 
expansions for the anomalous QCD dimensions 7 and 7^ defined as 

/y ry 

l(y')dy'), (n F )(xexp(J ^ F (y')dy') (16) 

are 

7 = 7o(l - ai7o - a 2 7o ~ a 37o) ( 17 ) 

with 7_p = 7 — r'/r, where 

r = (no)/ (n F ) = r (l - ri7 - r 2 ^l - r 3 ^). (18) 

All the coefficients a^Vi have been calculated and tabulated in J7|. The integrals Vj and 
terms in the right-hand sides of equations (fH]) , (115]) proportional to (n G ) 2 z 2 /2 are given in 
the Appendix. The corresponding terms in the left-hand sides can be written as 



[lnG G ]"= *^^[1 + E« r pnG F r = ^^1 + E«, (19) 



where 

M x = -(201 + 0.55 + 4/0, 

M 2 = aj- 2a 2 + 5 Gl - 4/ 2 + 4/i(2ai + 1.55), 

M 3 = 2 ai a 2 - 2a 3 + 0.55(3a 2 - 50 - 4/ 3 + 2/ 2 (4ai + 55) 

-4/! (a? - 2a 2 + 25a! + 0.755 2 ), 

3 

iVi = -2 fll - 0.55 - (1 + — )0i + 2ri, 
3 3 

iV 2 = a 2 - 2a 2 + 5ai - (1 + — )0 2 + (1 + — )0i(2ai + 1.55 - 2n) + 
r r 

2r 2 + n(3n -4ai -35), 
3 3 

iV 3 = 2 ai a 2 - 2a 3 + 0.55(3a 2 - 50 - (1 + — )fa + (1 + — )M 2a i + 

'n 'o 

2.55 - 2r0 - (1 + — )^i[a? - 2a 2 + 25 fll + 0.755 2 + 
ro 

2r 2 + n(3ri - 4ai - 55)] + 2r 3 + r 2 (6r! - 4a x - 55) + 
ri(4r 2 - 7.55ri - 6 ai n + 45 fll + 2a\ - 4a 2 + 1.55 2 ). (20) 

The terms with the same power of 70 in expressions on both sides should be equal. Therefore, 
one gets 



/l = i[4v 1 -2a 1 -B/2-^(l-A + «]] > (21) 



0i = 1 + 3 /ro (/i + 2r i + 2v 7 " 2fl i " 5 / 2 )> ( 22 ) 
/ 2 = -/i(4ai - 2vi + 35) + i(a 2 - 2a 2 - 4v 2 - 4a lVl + 5 Ql - 45v0 



n f \ ^ 1 A 1 \ 2 X ri ,1 1 , , m ,l 5 2 

, V 5/-, ^0 

3 

3„ _ , 1 



4-1(1 -v)^ (23) 



2 = 0i (2 fll + ^5 - 2n) + 1 + 3 /rQ (/2 + 2r 2 + 3r 2 - 2a 2 - 4am 

4 

+a\ - 35ri + 5 Ql - 2v 7 (f 1 + Ql + 5) + 4v 8 + — v 10 ) (24) 



h = oi d f + Y. d i Y i + rf 12 V 12 + 43V13], (25) 
6 i=l 

where df = M 3 + 4/ 3 , 

/ 3\-i 11 
3 = 1 + — [h + d^ + Y, d&i + d u v u ] , (26) 

where d^ = N 3 + (1 + ^)03- The coefficients dj are given in the Appendix. The numerical 
values of 0, for different number of active flavors are shown in the Table 1. 

Table 1 



n f 


ft 


h 


h 




02 


03 


3 


0.364 


-0.0279 


0.795 


0.637 


-0.276 


2.12 


4 


0.358 


-0.0457 


0.740 


0.631 


-0.286 


2.04 


5 


0.352 


-0.0629 


0.689 


0.625 


-0.295 


1.95 


S 


0.313 


0.310 


-0.120 


0.313 


0.310 


-0.120 



Herefrom it is easy to see that the asymptotical (70 — > 0) values of F G and are 
different as has been known since long ago ||10|| : 

F G - - F F - 1 4- — - - (271 

r 2,as — ^2,05 _ 3 _ 4' ^ ' 

The experimental values for 41.8 GeV gluon jets, F 2 G = 1.023 ± 0.008 ± 0.011, and for 
45.6 GeV uds quark jets, F% = 1.0820 ± 0.0006 ± 0.0046, are much lower than the above 
asymptotical limits. If one accepts, however, the effective value of as, averaged over all the 
energies of partons during the jet evolution, to be 0.2, then one gets F^(NLO) ~ 1.039 
and F^iNLO) ~ 1.068, by taking into account only the first correction proportional to 
70 . This is quite close to the experimental results 0. For as = 0.12, one gets about 10% 
higher values of i^'s. Thus we conclude that the NLO-approximation describes the widths 
at Z°-energies within 10-15% accuracy. 

At lower energies the widths should be slightly smaller due to the slow increase of as 
and somewhat smaller effective values of n/ leading to larger /1, 0i. Using the values of as 
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given in the PDG-data |TT[ and rif = 4, we predict the energy dependence of NLO- values of 
second factorial moments shown by the solid curves in Fig. 1. 

The additional dots at the ends of these curves demonstrate the effect due to possible change 
of effective values of rif to 3 at Q = 10 GeV and to 5 at Q = 90 Gev, where Q is the total 
energy as in |TT|. They show how small is the indefmiteness imposed by effective number of 



active flavours which is the only free parameter in such an approach . The curves indicate 
that the widths are closer to Poissonian ones at lower energies. Qualitatively, it agrees with 
experimental trends observed by the DELPHI collaboration 0. The coupling strength as 
changes within the interval of Q shown here from 0.18 to 0.12. The choice of the energy scale 
for jets is, however, highly nontrivial (see, e.g., @||). Therefore we do not plot experimental 
values here, just claiming the qualitative agreement within 15% accuracy. 

The cut-off of the integration region at e = e~~ y ~ e -^/Poa s f r0 m below and at 1 — e 
from above is not very crucial at present energies as seen from the dashed curves in Fig. 1. 
It diminishes the correction terms and, therefore, slightly increases the widths and flattens 
their energy dependence. Thus, the role of the power corrections is not very important. 

Unfortunately, the higher order corrections do not improve our estimates. On the con- 
trary, the 2NLO-term is positive and tends to violate slightly the agreement with experiment 
while 3NLO corrections are negative and so large that lead even to sub-Poissonian widths 
of distributions (F 2 < 1) for as = 0.2. These terms are approximately equal to NLO correc- 
tions due to large values of f 3 and 3 . The inception of such large values can be traced to 
rather large contributions of integrals containing In 2 x, i.e. to the region of very soft gluons. 
Thus the cut-off at e~ y and 1 — e~ y becomes more important for these terms. 

The 3NLO-corrections are overestimated due to the adopted Taylor series expansion with 
the assumption y ^> | lnx| which is invalid for soft gluons. For example, the k t - dependence 
of the coupling strength is transformed so that 

1 1 / lnx(l — x)\ , . 

y + la.x{l-x) y \ y J 

and the second term becomes infinitely large at x — > 0. The cut-off at x — e~ y leads to 



a factor of 2 only. Thus the above expansion implies some special presumption about the 
coupling strength behavior in the nonperturbative region as well as its modification at the 
limits of the perturbative one. The series (|T2|), (|i~3|) are with the sign-changing and increasing 
(in modulus) terms. From the values of fi,4>i in Table 1, one concludes that higher order 
terms are more important for the width of a quark jet compared with a gluon jet. 

The slopes of the widths are especially sensitive to these higher order terms because 
each of them is enlarged by the factor n when differentiating 7^. Thus 3NLO contribution 
is about 3 times larger than the NLO term in the slopes of widths. It demonstrates that 
any precise quantitative estimates of slopes become impossible. In particular, at present 
energies one cannot trust NLO estimates of these slopes as being small: Ftf'(NLO) ~ 
0.04; F 2 F '(NLO) ~ 0.092 at as = 0.2. However, one can predict the asymptotical value of 
the ratio of slopes as 

[Fl)' as ~ 210, ~ °' 43 ' (2J) 

which surely coincides with their NLO ratio. It demonstrates that the second factorial 
moment of quark jets approaches its asymptotical value faster than for gluon jets. 

Let us stress that all slopes and curvatures in pQCD are related to the running property 
of the QCD coupling constant since they are proportional to its derivatives which are equal 
to zero for a fixed coupling constant. It is interesting to note that the two- loop term in 
as proportional to Pi contributes only to the left-hand side of Eqs (|T9"D, namely, to the 
coefficients M3 and N3 in (|20"|) , and its role is very mild there (about 1-2 %). Thus it can 
be accounted with high accuracy considering it only in the expressions for the expansion 
parameter 70 . 

One can also easily check that all the relations of SUSY QCD (where nf = N c = Cp) are 
valid for all the coefficients shown above (e.g., F 2 G = etc.). The SUSY values of /j = 4>i 
are shown in the lower line of the Table 1 marked by S. The asymptotical SUSY values of 
F 2 are equal to 4/3. 
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III. THIRD MOMENTS OF THE MULTIPLICITY DISTRIBUTIONS 



The system of equations flnj) , (|T3) valid up to 3NL0-terms of 7q was also applied by us 
to calculation of third moments of the multiplicity distributions. It was done by equating 
the terms proportional to z 3 on both sides of the equations and writing down the third 
factorial moments (defined by Eq. (0) at q — 3) as 



3 3 

1=1 1=1 



(30) 



Proceeding in the same way as done above, one gets h Q = 9/4, g Q = 1 + r + Tq/4 and 
the values of the coefficients hi,gi in ( j30|) shown in the Table 2. The asymptotic limit of 
the third moment of quark jets is about twice larger than that of gluon jets. The analytic 
expressions are too lengthy to be presented here. 

Table 2 



n f 


hi 


h 2 


h 


9i 


92 


93 


3 


0.986 


-0.342 


2.49 


1.61 


-1.58 


7.74 


4 


0.972 


-0.380 


2.36 


1.60 


-1.59 


7.54 


5 


0.957 


-0.417 


2.25 


1.59 


-1.60 


7.34 


S 


0.844 


0.722 


-1.09 


0.844 


0.722 


-1.09 



Comparing f), <pi with /ij, gi one concludes that the corrections increase for higher mo- 
ments even in the NLO-approximation. Moreover, at present values of 70 ~ 0.5 they are 
rather large. Let us note the similarity in the structure of corrections for widths and third 
moments. They alternate in sign, and third coefficients are larger than the first ones. It is an 
indication on the sign-alternating asymptotic series, and Borel summation can be effective 
here. The increase of the coefficients originates from the terms containing the integrals of 
the type J In™ xdx oc n\. The termination of the cascade at e = e~ v leading to power cor- 
rections becomes more important below Z°. This reminds of the situation with renormalons 
(see, e.g., @). 
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In SUSY QCD the asymptotical values of third moments are equal to 9/4. The first 
correction given by /ii(SUSY)=g 1 (SUSY)=0.844 is almost as large as for ordinary gluon 
jets. It is similar to the correction for second moments. However, NNLO and 3NLO-terms 
for moments in SUSY QCD differ drastically from those for ordinary jets both in absolute 
values and signs as seen from the Tables 1 and 2. It demonstrates their sensitivity to the 
value of tq which is drastically different in the two cases. 

The similar procedure can be used for higher rank moments as well. 



IV. GENERAL DISCUSSION AND CONCLUSIONS 

The equations (1), (2) are dealing with probabilities and, therefore, are of classical nature. 
However, the quantum-mechanical interference has been accounted in the angular ordering 
effect. Nevertheless, there is no proof of their validity at all orders. The approach advocated 
above treats these equations as the kinetic equations in QCD for partonic cascade processes. 
Implicitly we have assumed that these equations describe the cascade down to extremely 
low energies of partons by imposing the limits and 1 of integration in the shares of energy. 
Thus the non-perturbative region of soft partons was assumed to be described in the same 
manner as the perturbative cascade. Probably, one should include into the consideration 
only the perturbative region by the requirement that the evolution parameters under the 
integral in (1), (2) stay always positive i.e. by introducing the cut-off at e and 1 — e 
where e = e~ y . Such a modification would lead to the power-like corrections which can be 
neglected asymptotically but contribute at present energies. Their role was considered for 
average multiplicities in 0. They are not very important for the moments in the NLO- 
approximation as shown above. However, more thorough treatment is needed, especially, in 
view of rather large contribution of soft gluons to f'3 and ^3. 

Other analytic solutions of these equations different from the perturbative one as well 



as nonperturbative modifications of the equations [IT3HTH] can be looked for. Especially 



interesting would be to learn more about the singularities of the generating functions in the 
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2-plane which are known up to now for the leading order solution only. 

Leaving this program for future studies, we can now compare the results obtained from 
perturbative solutions of the equations for average multiplicities [3-7] and for widths of the 
multiplicity distributions of gluon and quark jets up to 3NLO approximation of pQCD. 
Leading order predictions for any quantity are quite far from present experimental data. 
NLO corrections are always pointing in the right direction of closer agreement with experi- 
ment. In particular, the energy dependence of the gluon jet average multiplicity, of the ratio 
of its slopes for gluon and quark jets and the values of widths at Z°-energies can be fitted 
with rather high accuracy. However, it still fails in the ratio of average multiplicities of gluon 
and quark jets differing from experiment by about 30% in absolute values and even more 
in its energy dependence slope 0. Here, 2NLO and 3NLO terms improve the situation so 
that the ratio r differs by about 10-15% only. They do not spoil good qualitative features 
of NLO in the energy dependence ||. However, these corrections become very large for the 
widths and for the third moments as shown above, as well as for the slope and curvature in 
the energy dependence of average multiplicities 0. 

Thus we conclude that present perturbative QCD results can describe the experimental 
data within 10-15% accuracy. The perturbative series breaks down, however, at different 
orders for different quantities. There seems to be no standard way to improve the analytic 
results and unanimously predict where one should truncate the expansion. Nevertheless, the 
general trends obtained from the perturbative approach are steadily indicating qualitative 
convergence of theory and experiment. Moreover, the computer results show that the 
exact solutions of QCD equations can be even closer to experiment. 
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APPENDIX: 



The terms proportional to (n G ) 2 z 2 /2 in the right-hand sides of equations flljj) , (|T^) up 
to 0(7q) corrections can be written as 

G G - 1 - if = |(1 - f ll0 - Mo 2 - /s7o 3 ), (Al) 
G G -> 2 7 if + ff = ^[1 - ( fll + A) 7 o + 

^(Oi/i- 02-/2 + 0.55/0], (A2) 

G G - 2(2 7 2 + 7 ')if + 4 7 if ' + if" = ^ 7o 2 [l - 7o (2 ai + /i + 0.55)], (A3) 

- 8 7 3 if = f 7o 3 , (A4) 

(g)^4 7o 3, (A5) 
(/i + ax) 



(S) / -^-5,o + 6)+ 7o2[ 



3 



4^?^(0i + a a - 2n) - 4 Ql (l - ro" 1 ) 2 - ^(1 - r 1 )} + 

2 7o 3 r - 2 [2(r - l)(2a 1 r 1 + 5r x - 2(r 2 + r 2 )) - 2a 2 (r - l) 2 + 2r 2 - 

4r 2 2 

_£ (ai/l _ fl2 _ /2 ) + + ro/3)(01 _ 2ri ) _ _/ ir g) + 

2(1 + r /3)(2r 2 - a 2 - 2 + 3r 2 + - 2am - tym)], (A6) 

(^)' - 8 7o 2 r - 2 (l - r /3) + 4 7o 3 r Q - 2 [4r 1 - ^ - 

(1 + r o /3)0! - (1 - r /3)(4oi + B)], (A7) 
f GfG'p\' 3 4(9-r ) 

1"gtJ 70 (A8) 



M ^7o 3 -' (A9) 

~* 7 ° 2 3r| + 7 ° 3 3^ [5n ~ 8ai ~ W ~ (1 + (A10) 

f^)'-7o 31 , (AH) 
\ G F J r 



GaG 



Gi 



7 3 4r - 2 (l + 5r /3), (A12) 



G 

G~g ' 3rg vu " u ' ~' ' ' u 3r 



- 1 - ^(1 ) - 5r o + 6) + 707r^(2/ir 2 - (3 + r o )0i + 4n(3 - r )), (A13) 
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GfG'p 

a 



G 



7o 4r - 2 (l-r /3), 
GqG' f 8 



Gi 



7o 



3r n 



The coefficients d,- in (|25T), (|26|) are 



d 1 = 2(2ai/ 1 -2o 2 -2/ 2 + 3S/ 1 ), 



d 2 = 2(4a x + 2.A + 55), 



d 3 = -4, 



" ? [r 2 / 2 - 0.5(3 + r o )0 2 - rg/ x (ai + 1.55) + 



2iV c r 2 



0.25(3 + r o )0i(2ai + 35 - 4r : ) + r 2 a 2 - 0.5(3 + r )a 2 + 2r 2 (3 - r ) + 
ri[r (2ai - 1.5ri + 35) - 6a x + 9n - 95] - 1.5a 2 (r - l) 2 ], 



27V c r 2 



[(3 + r o )0i + 4(3 - r )ai + 5(21 - 13r + 2r 2 ) - n(12 - r )], 



dp, 



n f (9 - r ) 
4iV c r 2 ' 



d 7 = 0.5di, 



dg — — d 2 , 



do = -d 3 



dio = -^[8^ + (1 + Sr^ 1 )^ + 25(3 + 2r ) - 5n] 



d 



ii 



3 + 5r 
2r n 2 ' 



13 



= -3, 



d 



3rii 



13 



4iV c r 2 



2 ' 



di4 = 3/r . 



The integrals Vj are as follows 

r 1 r 1 3 

vi = / Vdx = / (1 — -x + x 
jo jo 2 

v 2 = / 
JO 



ln(l — x) 



2T/lnx(l -x) 



x\ , 11 
— lax = — 
2 ; 24 

67 - 6tt 2 



dx = 



36 



v 3 



= /' 

JO 



ln 2 (l - x) 

x 



2\/(ln 2 x + ln 2 (l -x)) 



dx = 2C(3) 



413 
108 



/•i 2 
v 4 = / [ x 2 + (1 - x) 2 ]rfx = - 

JO 

v 5 = / [x 2 + (1 — x) 2 ] lnxcix = 

JO 

v 6 = J [x 2 + (1 — x) 2 ] In 2 xdx 

/(i-f)«H 



v 7 = / $dx 
/o 



v 8 

V 9 : 



13 

18 
89 

54 
3 



JO 



1 7 
$ In xdx = — 

8 

$ In 2 xdx = — 



Vio 



Vn 



L 



i 

L 



i r 

L 



• -i 

X 

1 



7T 



2 5 



x 



ln(l — x)dx = 

V ' 6 8 

ln 2 (l-x)rfx = |-2C(3) 



V12 = / (x _1 - 210 lnxln(l - x)rfx = C(3) - 395/216 + 11tt 2 /72 
jo 

Vl3 = / [x 2 + (l-x) 2 ]lnxln(l-x)rfx = 71/54 -tt 2 /9 

JO 

vu= / (l-0.5x-x- 1 )lnxln(l-x)rfx = 1.5-tt 2 /8-C(3) 

JO 



(A27) 
(A28) 
(A29) 



(A30) 
(A31) 

(A32) 

(A33) 
(A34) 
(A35) 
(A36) 
(A37) 
(A38) 
(A39) 
(A40) 
(A41) 
(A42) 
(A43) 



£ means Riemann's ^-function. 
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FIG. 1. The energy behaviour of the second factorial moments of quark (F) and gluon (G) jets. 
The limits of integration are chosen as and 1 (solid lines) or e and 1 — e (dashed lines). The dots 
at the ends show that the curves are insensitive to variation of the effective number of flavors (see 
text). 
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